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ABSTRACT 

We investigate the properties of relativistic r-modes of slowly rotating neu- 
tron stars by using a relativistic version of the Cowling approximation. In our 
formalism, we take into account the influence of the Coriolis like force on the 
stellar oscillations, but ignore the effects of the centrifugal like force. For three 
neutron star models, we calculated the fundamental r-modes with V = m = 2 
and 3. We found that the oscillation frequency a of the fundamental r-mode is 
in a good approximation given by a ~ k,q Q, where a is defined in the corotating 
frame at the spatial infinity, and Q is the angular frequency of rotation of the star. 
The proportional coefficient kq is only weakly dependent on Q, but it strongly 
depends on the relativistic parameter GM/c 2 R, where M and R are the mass 
and the radius of the star. All the fundamental r-modes with V = m computed 
in this study are discrete modes with distinct regular eigenfunctions, and they all 
fall in the continuous part of the frequency spectrum associated with Kojima's 
equation (Kojima 1998). These relativistic r-modes are obtained by including 
the effects of rotation higher than the first order of Q so that the buoyant force 
plays a role, the situation of which is quite similar to that for the Newtonian 
r-modes. 

Subject headings: instabilities — stars: neutron — stars: oscillations — stars: 
rotation 

1. Introduction 

It is Andersson (1998) and Friedman & Morsink (1998) who realized that the r-modes 
in rotating stars are unstable against the gravitational radiation reaction. Since then, a large 
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number of papers have been published to explore the possible importance of the instability 
in neutron stars. The r-modes in rotating stars are restored by the Coriolis force. They have 
the dominant toroidal component of the displacement vector and their oscillation frequencies 
are comparable to the rotation frequency Q of the star (see, e.g., Bryan 1889; Papaloizou 
k Pringle 1978; Provost, Berthomieu, & Roca 1981; Saio 1982; Unno et al. 1989). For the 
r-mode instability, see reviews by, e.g., Friedman & Lockitch (1999), Andersson & Kokkotas 
(2001), Lindblom (2001), and Friedman & Lockitch (2001). Most of the studies on the r-mode 
instability in neutron stars, however, have been done within the framework of Newtonian 
dynamics. Since the relativistic factor can be as large as GM/c 2 R ~ 0.2 for neutron stars 
where M and R are respectively the mass and the radius, the relativistic effects on the 
r-modes are essential. 

From time to time the effects of general relativity on the r-modes in neutron stars have 
been discussed in the context of the r-mode instability. Kojima (1998) is the first who 
investigated the r-modes in neutron stars within the framework of general relativity. In 
the slow rotation approximation, Kojima (1998) derived a second order ordinary differential 
equation governing the relativistic r-modes, expanding the linearized Einstein equation to 
the first order of Q, and assuming that the toroidal component of the displacement vector is 
dominant and the oscillation frequency is comparable to Q. Kojima (1998) showed that this 
equation has a singular property and allows a continuous part in the frequency spectrum of 
the r-modes (see, also, Beyer & Kokkotas 1998). Recently, Lockitch, Andersson, & Friedman 
(2001) showed that Kojima's equation is appropriate for non-barotropic stars but not for 
barotropic ones. They found that both discrete regular r-modes and continuous singular 
r-modes are allowed in Kojima's equation for uniform density stars, and suggested that the 
discrete regular r-modes are a relativistic counterpart of the Newtonian r-modes. Yoshida 
(2001) and Ruoff & Kokkotas (2001a) showed that discrete regular r-mode solutions to 
Kojima's equation exist only for some restricted ranges of the polytropic index and the 
relativistic factor for polytropic models, and that regular r-mode solutions do not exist for 
the typical ranges of the parameters appropriate for neutron stars. 

The appearance of continuous singular r-mode solutions in Kojima's equation (Kojima 
1998) caused a stir in the community of people who are interested in the r-mode instability. 
Lockitch et al. (2001) (see also Beyer & Kokkotas 1999) suggested that the singular property 
in Kojima's equation could be avoided if the energy dissipation associated with the gravita- 
tional radiation is properly included in the eigenvalue problem because the eigenfrequencies 
become complex due to the dissipations and the singular point can be detoured when inte- 
grated along the real axis. Very recently, this possibility of avoiding the singular property in 
Kojima's equation has been examined by Yoshida & Futamase (2001) and Ruoff & Kokkotas 
(2001b), who showed that the basic properties of Kojima's equation do not change even if the 
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gravitational radiation reaction effects are approximately included into the equation. Quite 
interestingly, as shown by Kojima & Hosonuma (2000), if the third order rotational effects 
are added to the original Kojima's equation (Kojima 1998), the equation for the r-modes 
becomes a fourth order ordinary linear differential equation, which has no singular properties 
if the Schwarzschild discriminant associated with the buoyant force does not vanish inside 
the star. By solving a simplified version of the extended Kojima's equation for a simple 
toy model, Lockitch & Andersson (2001) showed that because of the higher order rotational 
terms the singular solution in the original Kojima's equation can be avoided. Obviously, we 
need to solve the complete version of the extended Kojima's equation to obtain a definite 
conclusion concerning the continuous singular r-mode solutions. 

It may be instructive to turn our attention to a difference in mathematical property 
between the Newtonian r-modes and the relativistic r-modes associated with Kojima's equa- 
tion. It is usually assumed that in the lowest order of Q the eigenfunction as well as the 
eigenfrequency of the r-modes is proportional to Q. In the case of the Newtonian r-modes, if 
we employ a perturbative method for r-modes in which the angular frequency Q is regarded 
as a small expanding parameter, the radial eigenfunctions of order of Q can be determined 
by solving a differential equation derived from the terms of order of Q 3 , which bring about 
the couplings between the oscillations and the buoyant force in the interior (e.g., Provost et 
al 1981, Saio 1982). In other words, there is no differential equation, of order of f2, which 
determines the radial eigenfunction of the Newtonian r-modes. On the other hand, in the 
case of the general relativistic r-modes derived from Kojima's equation, we do not have to 
take account of the rotational effects of order of fi 3 to obtain the radial eigenfunctions of 
order of Q. That is, the eigenfrequency and eigenfunction of the r-modes are both deter- 
mined by a differential equation (i.e., Kojima's equation) derived from the terms of order of 
Q. This remains true even if we take the Newtonian limit of Kojima's equation to calculate 
the r-modes (Lockitch et al. 2001; Yoshida 2001). We think that this is an essential differ- 
ence between the Newtonian r-modes and the relativistic r-modes associated with Kojima's 
equation. Considering that Kojima's equation can give no relativistic counterpart of the 
Newtonian r-mode, it is tempting to assume that some terms representing certain physical 
processes are missing in the original Kojima' equation. On the analogy of the r-modes in 
Newtonian dynamics, we think that the buoyant force in the interior plays an essential role 
to obtain a relativistic counterpart of the Newtonian r-modes and that the terms due to 
the buoyant force will appear when the rotational effects higher than the first order of Q 
are included. This is consistent with the suggestions made by Kojima & Hosonuma (2000) 
and Lockitch & Andersson (2001). In this paper, we calculate relativistic r-modes by taking 
account of the effects of the buoyant force in a relativistic version of the Cowling approxi- 
mation, in which all the metric perturbations are omitted. In our formulation, all the terms 
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associated with the Coriolis force are included but the terms from the centrifugal force are all 
ignored. Our formulation can take account of the rotational contributions, due to the Coriolis 
force, higher than the first order of Q. Note that our method of solution is not a perturba- 
tion theory in which Q is regarded as a small expanding parameter for the eigenfrequency 
and eigenfunction. Similar treatment has been employed in Newtonian stellar pulsations in 
rotating stars (see, e.g., Lee & Saio 1986; Unno et al. 1989; Bildsten, Ushomirsky & Cutler 
1996; Yoshida & Lee 2001). This treatment is justified for low frequency modes because the 
Coriolis force dominates the centrifugal force in the equations of motion. The plan of this 
paper is as follows. In §2, we describe the formulation of relativistic stellar pulsations in the 
relativistic Cowling approximation. In §3, we show the modal properties of the r-modes in 
neutron star models. In §4, we discuss about our results, specially the effect of the buoyancy 
in the fluid core on the modes. §5 is devoted to conclusions. In this paper, we use units in 
which c — G—1, where c and G denote the velocity of light and the gravitational constant, 
respectively. 



2. Formulation 

2.1. Equilibrium State 

We consider slowly and uniformly rotating relativistic stars in equilibrium. If we take 
account of the rotational effects up to first order of fi, the geometry in the stars can be 
described by the following line element (see, e.g. Thorne 1971): 

ds 2 = g aP dx a dx p = -e 2v{r) dt 2 + e 2X{r) dr 2 + r 2 d6 2 + r 2 sin 2 6d(p 2 - 2u{r) r 2 sin 2 Odtdip . (1) 

The fluid four-velocity in a rotating star is given by 

u a = j(r,9)(t a + n V a ), (2) 

where t a and </? Q stand for the timelike and rotational Killing vectors, respectively. Here the 
function 7 is chosen to satisfy the normalization condition u a u a — — 1. If we consider the 
accuracy up to order of Q, the function 7 reduces to: 

7 = e- u{r) . (3) 

Once the physical quantities of the star such as the pressure, p(r), the mass-energy density, 
p(r), and the metric function, v(r), are given, the rotational effect on the metric, u(r) can 
be obtained from a well-known numerical procedure (see, e.g., Thorne 1971). 
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2.2. Pulsation Equations in the Cowling Approximation 

General relativistic pulsation equations are usually obtained by linearizing Einstein's 
field equation. The linearized Einstein equation contains perturbations associated with the 
metric fluctuations and the fluid motions. In this paper, to simplify the problem, we employ 
a relativistic version of the Cowling approximation, in which all the metric perturbations 
are omitted in the pulsation equations (see McDermott, Van Horn, & Scholl 1983, and Finn 
1988). The relativistic Cowling approximation is accurate enough for the /- and p-modes 
in non-rotating stars (Lindblom & Splinter 1992). It is also the case for the modes in 
slowly rotating stars (Yoshida & Kojima 1997). The relativistic Cowling approximation is 
a good approximation for oscillation modes in which the fluid motions are dominating over 
the metric fluctuations to determine the oscillation frequency. Therefore, it is justified to 
employ the relativistic Cowling approximation for the r-modes, for which the fluid motion 
is dominating. 

If we employ the Cowling approximation, we can obtain our basic equations for pulsa- 
tions from the perturbed energy and momentum conservation laws: 

5 {u a V pT^ ) = , (energy conservation law) (4) 
5 (g"V pT^) = , (momentum conservation law) (5) 

where V a is the covariant derivative associated with the metric, Tg is the energy-momentum 
tensor, and q^ is the projection tensor with respect to the fluid four-velocity. Here, 5Q 
denotes the Eulerian change in the physical quantity Q. In this paper, we adapt the adiabatic 
condition for the pulsation: 

Ap=^—Ap, (6) 
P + P 

where T is the adiabatic index defined as 

r = e ± E (%) , (7) 

ad 



and AQ stands for the Lagrangian change in the physical quantity Q. The relation between 
the Lagrangian and the Eulerian changes is given by the equation: 

AQ = 5Q + £ C Q , (8) 

where £^ is the Lie derivative along the Lagrangian displacement vector ( a , which is defined 
by the relation: 

5u a = q a ^ = q^£ u Cf- (9) 
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Notice that we have 5u a = 5u a in the Cowling approximation. Because we are interested 
in pulsations of stationary rotating stars, we can assume that all the perturbed quantities 
have time and azimuthal dependence given by e tcrt+tmip , where m is a constant integer and 
o is a constant frequency measured by an inertial observer at the spatial infinity. Because 
of this assumption, the relation between the Lagrangian displacement ( a and the velocity 
perturbation 5u a reduces to an algebraic relationship: 

5u a = i 1 a( a , (10) 

where a is the frequency defined in the corotating frame defined as a = a + mfl. Note that 
the gauge freedom in ( a has been used to demand the relation u a ( a = 0. By substituting 
equations (6) and (10) into equations (4) and (5), we can obtain the perturbed energy 
equation, 

- v a ( 7 C a ) + ^(Sp + C Q v a p) = o , (ii) 

and the perturbed momentum equation, 

-^a 2 g a ^ + 2i 7 aC^ a + <£S7 P (-^-) + (-^- € + ^f) A, = 0, (12) 

\P + Pj \P + P P + P/ 

where A a is the relativistic Schwarzschild discriminant defined by 

A a = V aP -^V a p. (13) 
p + p Tp 

Notice that equations (11) and (12) have been derived without the assumption of slow 
rotation. Physically acceptable solutions of equations (11) and (12) must satisfy boundary 
conditions at the center and the surface of the star. The surface boundary condition at 
r = R is given by 

Ap = 5p + CV a p = 0, (14) 

and the inner boundary condition is that all the eigenfunctions are regular at the center 
(r = 0). 

On the analogy between general relativity and Newtonian gravity, the second term 
on the left hand side of equation (12) is interpreted as a relativistic counterpart of the 
Coriolis force. In our formulation, the terms due to the Coriolis like force are included in the 
perturbation equations, but the terms due to the centrifugal like force, which are proportional 
to Vl 2 /(GM/R 3 ), are all ignored. In Newtonian theory of oscillations, this approximation is 
justified for low frequency modes satisfying the conditions \2Q/a\ > 1 and VL 2 / (GM/R 3 ) <^ 1 
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(Lee & Saio 1986; Unno et al. 1989; Bildsten, Ushomirsky & Cutler 1996; Yoshida & Lee 
2001). In general relativity, we note that it is difficult to make a clear distinction between 
inertial forces such as the Coriolis force and the centrifugal force. From a physical point of 
view, however, it is also acceptable to use the approximation for low frequency oscillations 
in the background spacetime described by the metric (1). 

The eigenfunctions are expanded in terms of spherical harmonic functions Y™^, <p) 
with different values of I for a given m. The Lagrangian displacement, ( k and the pressure 
perturbation, 5p/ (p + p) are expanded as 

oo 

C = rJ2 Si(r)YT (0, <p)e iat , (15) 

l>\m\ 



oo 



J,/'>|m| k r ' 

^ £ W^^r)^™}.", (17 ) 



sin 2 # ^ , l W fy> w <9fl 



5p 



E^MWrf^, (18) 

i>|m| 

where / = |m| + 2k and /' = I + 1 for even modes and / = \m\ + 2k + 1 and /' = / — 1 for odd 
modes where A; = 0, 1, 2 • • • (Regge & Wheeler 1957; Thorne 1980). Here, even and odd 
modes are, respectively, characterized by their symmetry and antisymmetry of the eigenfunc- 
tion with respect to the equatorial plane. Substituting the perturbed quantities (15)-(18) 
into linearized equations (11) and (12), we obtain an infinite system of coupled ordinary 
differential equations for the expanded coefficients. The details of our basic equations are 
given in the Appendix. Note that non-linear terms of q = 2u/a where Co = Q — uj are kept 
in our basic equations. For numerical calculations, the infinite set of ordinary differential 
equations are truncated to be a finite set by discarding all the expanding coefficients asso- 
ciated with I larger than Z max , the value of which is determined so that the eigenfrequency 
and the eigenfunctions are well converged as / max increases (Yoshida & Lee 2000a). 



3. r-Modes of Neutron Star Models 

The neutron star models that we use in this paper are the same as those used in the 
modal analysis by McDermott, Van Horn, & Hansen (1988). The models are taken from 
the evolutionary sequences for cooling neutron stars calculated by Richardson et al. (1982), 



- 8- 



where the envelope structure is constructed by following Gudmundsson, Pethick & Epstein 
(1983). These models are composed of a fluid core, a solid crust and a surface fluid ocean, 
and the interior temperature is finite and is not constant as a function of the radial distance 
r. The models are not barotropic and the Schwarzschild discriminant \A\ has finite values in 
the interior of the star. In order to avoid the complexity in the modal properties of relativistic 
r-modes brought about by the existence of the solid crust in the models (see Yoshida & Lee 
2001), we treat the whole interior of the models as a fluid in the following modal analysis. 

We computed frequency spectra of r-modes for the neutron star models called NS05T7, 
NS05T8, and NS13T8 (see, McDermott et al. 1988). The physical properties such as the 
total mass M, the radius R, the central density p c , the central temperature T c and the 
relativistic factor GM/c 2 R are summarized in Table 1 (for other quantities, see McDermott 
et al. 1988). In Figures 1 and 2, scaled eigenfrequencies k = a/tt of the r-modes of the 
three neutron star models are given as functions of Q = Q/ ^/GM/R 3 for m = 2 and 3 cases, 
respectively. Here only the fundamental r-modes with V — m are considered because they 
are most important for the r-mode instability of neutron stars. We note that it is practically 
impossible to correctly calculate rotationally induced modes at Cl ~ because of their 
coupling with high overtone g-modes having extremely low frequencies. From these figures, 
we can see that the scaled eigenfrequency k is almost constant as tl varies. In other words, 
the relation a ~ k q VL is a good approximation for the fundamental r-modes with V = m, 
where kq is a constant. Comparing the two frequency curves, which nearly overlap each 
other, for the models NS05T7 and NS05T8, it is found that the detailed interior structure 
of the stars such as the temperature distribution T(r) does not strongly affect the frequency 
of the fundamental r-modes with /' = m. This modal property is the same as that found 
for the fundamental V — m r-modes in Newtonian dynamics (see, Yoshida & Lee 2000b). 
On the other hand, comparing the frequency curves for the models NS05T7 (NS05T8) and 
NS13T8, we note that the r-mode frequency of relativistic stars is strongly dependent on the 
relativistic factor GM/c 2 R of the models. This is because the values of the effective rotation 
frequency uj = Q — u in the interior is strongly influenced by the relativistic factor. Similar 
behavior of the GM/c 2 R dependence of the r-mode frequency has been found in the analysis 
of Kojima's equation (Yoshida 2001; Ruoff & Kokkotas 2001a). In Table 2, the values of k 
for the r-modes shown in Figures 1 and 2 are tabulated, where k,q are evaluated at Cl — 0.1. 
The boundary values for the continuous part of the frequency spectrum for the /' = m = 2 
r-modes derived from Kojima's equation (Kojima 1998) are also listed in the same table. 
As shown by Kojima (1998) (see, also, Beyer & Kokkotas 1999; Lockitch et al. 2001), if an 
r-mode falls in the frequency region bounded by the boundary values, Kojima's equation 
becomes singular and yields a continuous frequency spectrum and singular eigenfunctions 
as solutions. Although all the r-modes obtained in the present study are in the bounded 
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frequency region, our numerical procedure shows that the r-modes obtained here are isolated 
and discrete eigenmodes. In fact, no sign of continuous frequency spectrum appears in the 
present numerical analysis (for a sign of the appearance of a continuous frequency spectrum 
in a numerical analysis, see Schutz & Verdaguer 1983). 

In Figures 3 and 4, the eigenfunctions % T 2 for the I' — m — 2 fundamental r-modes in 
the neutron star models NS05T8 and NS13T8 at Cl — 0.1 are shown. We can confirm from 
these figures that the eigenfunctions show no singular property, even though the frequencies 
are in the continuous part of the spectrum associated with Kojima's equation. These figures 
also show that the fluid motion due to the r-modes is more confined near the stellar surface 
for NS13T8 than for NS05T8. This suggests that the eigenfunctions iT m tend to be confined 
to the stellar surface as the relativistic factor of the star increases. The same property 
appears in regular r-mode solutions derived from Kojima's equation (Yoshida 2001; Ruoff 
& Kokkotas 2001a; Yoshida & Futamase 2001). But, in the present case, the confinement 
of the eigenfunction for NS13T8 is not so strong, even though the model NS13T8 is highly 
relativistic in the sense that the relativistic factor is as large as GM/c 2 R = 0.249. 



4. Discussion 

In the neutron star models analyzed in the last section, the buoyant force is produced by 
thermal stratification in the star. However, as suggested by Reisenegger & Goldreich (1992), 
the buoyant force in the core of neutron stars might become even stronger if the effect of the 
smooth change of the chemical composition of charged particles (protons and electrons) in the 
core is taken into account. We thus think it legitimate to examine the effects of the enhanced 
buoyant force on the modal properties of the relativistic r-modes. Although our neutron 
star models do not provide enough information regarding the composition gradient, we may 
employ for this experiment an approximation formula for the Schwarzschild discriminant due 
to the composition gradient in the core given by Reisenegger & Goldreich (1992): 

rA r = 3.0 x 10- 3 (19) 
\pnucj par 

where p nuc denotes the nuclear density p nuc = 2.8 x 10 14 gcm~ 3 . In Figure 5, Brunt- Vaisala 
frequency due to the thermal stratification (solid line) and that due to the composition 
gradient (dashed line) are given as a function of log(l — r/R) for model NS13T8, where 
the frequencies are normalized by (GM/i? 3 ) 1 / 2 . In this paper, the relativistic Brunt- Vaisala 
frequency N is defined as 

N 2 = -A_ $ . (20) 
p + p dr 
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This figure shows that in the core the Brunt- Vaisala frequency due to the composition 
gradient is by several orders of magnitude larger and hence the corresponding buoyant force 
is much stronger than those produced by the thermal stratification. Note that the Brunt- 
Vaisala frequency due to the composition gradient becomes as large as ~ (GM/i? 3 ) 1 / 2 in the 
core (see also, e.g., Lee 1995). 

Replacing the Schwarzschild discriminant in the original neutron star models with that 
due to the composition gradient calculated by using equation (19), we computed the funda- 
mental r-modes with I' = m = 2 for the models NS05T8 and NS13T8, and plotted k = cr/Cl 
against rotation frequency fl in Figure 6, where «'s for the original models with the thermal 
stratification were also plotted for convenience. As shown by the figure, k slightly increases 
as ti — > for the models with the enhanced buoyant force. This is remarkable because 
such behavior of k for the fundamental r-modes with I' — m is not found in the case of 
Newtonian r-modes in non-barotropic stars (Yoshida & Lee 2000b). The slight increase of 
k with decreasing fl may be explained in terms of the property of the eigenfunctions of the 
r-modes. We depict the eigenfunctions 1T2 of the V — m — 2 fundamental r-modes in the 
neutron star model NS13T8 at Cl — 0.2 and .02 in Figures 7 and 8, where the solid line 
and dashed line in each figure denote the r-modes in the models with the compositional 
stratification and with the thermal stratification, respectively. As shown by Figure 7, in the 
case of rapid rotation, the buoyant force in the core does not strongly affect the properties 
of the eigenfunctions. On the other hand, in the case of slow rotation, Figure 8 shows that 
the amplitude of the eigenfunction is strongly confined to the region near the surface for the 
model with the enhanced buoyant force. For relativistic r-modes, as a result of the general 
relativistic frame dragging effect, the effective rotation frequency uo = Q — uo oc Q acting on 
local fluid elements is an increasing function of the distance r from the stellar center, and 
thus fluid elements near the stellar surface are rotating faster than those near the stellar 
center. The oscillation frequency of the r-modes may be determined by the mean rotation 
frequency obtained by averaging local rotation frequencies over the whole interior of the 
star with a certain weighting function associated with the eigenfunction. In this case, it is 
reasonable to expect that the relativistic r-modes that have large amplitudes only in the 
regions near the surface get larger values of k than those that have large amplitudes deep in 
the core. 

For the models with the enhanced buoyant force, the current quadrupole moment of the 
r modes that determines the strength of the instability will be largely reduced when fl is small 
as a result of the strong confinement of the amplitudes. But, since the r-mode instability is 
believed to operate at rapid rotation rates, the effect of the amplitude confinement at small 
rotation rates may be irrelevant to the instability. 
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5. Conclusion 

In this paper, we have investigated the properties of relativistic r-modes in slowly ro- 
tating neutron stars in the relativistic Cowling approximation by taking account of higher 
order effects of rotation than the first oder of Q. In our formalism, only the influence of the 
Coriolis like force on the oscillations are taken into account, and no effects of the centrifugal 
like force are considered. We obtain the fundamental r-modes associated with I' — m — 2 
and 3 for three neutron star models. We find that the fundamental r-mode frequencies 
are in a good approximation given by a ~ k,q ^- The proportional coefficient k,q is only 
weakly dependent on Q, but strongly depends on the relativistic parameter GM/c 2 R. For 
the fundamental r-modes with V = m, we find that the buoyant force in the core is more 
influential to the relativistic r-modes than to the Newtonian r-modes. All the r-modes ob- 
tained in this paper are discrete modes with distinct regular eigenfunctions, and they all fall 
in the frequency range of the continuous spectrum of Kojima's equation. We may conclude 
that the fundamental r-modes obtained in this paper are the relativistic counterpart of the 
Newtonian r-modes. 

Here, it is legitimate to mention the relation between the present work and the study 
by Kojima & Hosonuma (1999), who studied relativistic r-modes in slowly rotating stars in 
the Cowling approximation. By applying the Laplace transformation to linearized equations 
derived for the r-modes, Kojima & Hosonuma (1999) examined how a single component of 
initial perturbations with axial parity evolves as time goes, and showed that the perturba- 
tions cannot oscillate with a single frequency. In this paper, we have presented a formulation 
for small amplitude relativistic oscillations in rotating stars in the relativistic Cowling ap- 
proximation, and solved the oscillation equations as a boundary-eigenvalue problem for the 
r-modes. In our formulation, we assume neither axially dominant eigenfunctions nor low 
frequencies to calculate the r-modes. Considering these differences in the treatment of the 
r-mode oscillations between the two studies, it is not surprising that our results are not 
necessarily consistent with those by Kojima & Hosonuma (1999). 

Our results suggest that the appearance of singular r-mode solutions can be avoided by 
extending the original Kojima's equation so that terms due to the buoyant force in the stellar 
interior are included. As discussed recently by Yoshida & Futamase (2001) and Lockitch & 
Andersson (2001), we believe that the answer to the question whether r-mode oscillations 
in uniformly rotating relativistic stars show true singular behavior may be given by solving 
the forth order ordinary differential equation derived by Kojima & Hosonuma (2000) for 
r-modes. Verification of this possibility remains as a future study. 

We would like to thank H. Saio for useful comments. We are grateful to the anonymous 
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A. Basic Equations 

We introduce column vectors yi, y 2 , h, and t, whose components are defined by 

yi,k = Si(r) , (Al) 
y2,k = ^;SU l (r), (A2) 

dr 

h, k = H l (r), (A3) 

and 

t, k = Mr), (A4) 

where / = \m\ + 2k — 2 and V — I + 1 for "even" modes, and I — \m\ + 2k — 1 and I' — I — 1 
for "odd" modes, and k — 1, 2, 3, 

In vector notation, equations for the adiabatic nonradial pulsation in a slowly rotating 
star are written as follows: 

The perturbed energy equation (11) reduces to 

r ^d7 + ( 3 ~ T + yi + T y2 ~ A ° h + ° 2 ^ & ^~ mh + = ' ^ 

The r component of the perturbed momentum equations (12) reduces to 

^ - (e 2A Cl a 2 + rA r ) Yl + (U + rA r )y 2 - Cl ff 2 x(-mh + C it) = . (A6) 
dr 



Here, 



u _r_iir£> y_J_hv (A7) 

£ ' dlnr' ( j 



e 2iy d 

x = r\ r " e 



r dr 

r 2 e~ 2v M , ,, ...">/ 

r 



2 , (A8) 

(7 / 



dr 

and (5 = lo/(GM/R 3 )^ 2 , a = a/(GM/R 3 ) 1/2 , and <r = a / (GM / R 3 ) 1 / 2 are frequencies in the 
unit of the Kepler frequency at the stellar surface, where cu = Q — u and a = a + mfl. 
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The 6 and ip components of the perturbed momentum equations (12) reduce to 
L h + M 1 it = ^2y 2 + {O + M 1 Lj- 1 K} fxyi + ^y 2 ) , 



where 



Lxit + M h = K ( xyi + ^^y 2 

V Ci<7 

O = mAo 1 - MiL^K . 



(A10) 
(All) 
(A12) 



The quantities C , K, L , L 1; A , A 1; M , Mi are matrices written as follows: 
For even modes, 

(Q))i,i = —(I + 2) J™i , (C )i+i,i = (I + 1)^+2 > 



(K) w = 
(Lo)i,i = 1 



J l+l 

Z + l' 
mq 



1(1 + 1)' 



(K)i ii+ i — 
(Li) M = 1 - 



Tin 
J l+2 

Z + 2 ' 



mq 



(/ + !)(/ + 2)' 



(A ) M = Z(Z + 1) , (Ai) M = (Z + 1)(Z + 2) , 



Z + 3 



(Mo) M = g— J^i , (M ),, )+1 = q— J™ 2 



Z + 2 



Z + l 



(Mi)y = ^ , ( M iW = ^ ^2 
where Z = |m| + 2i — 2 for « = 1, 2, 3, . . ., and q = 2u/a, and 



J, m = 



(Z + m)(l — m) 
(21 - 1)(2Z + 1) 



1/2 



For odd modes, 



(Q))i,i — (Z — l)^ m 5 (Co)i,i+i — — (Z + 2) J™ x , 



(Lo) 



jm 

(K) M = — — , 
mg 



;k) 



H-l 



Z + l' 



mg 



Z(Z + 1)' 
(Ao) M = Z(Z + 1) , 

(Mo),, = g^ ^ , 



(Li)m " 1 " 1(1 -iy 

(Ai) M = Z(Z - 1) , 



(A13) 



-14- 



where / = |m| + 2i — 1 for i — 1, 2, 3, . . .. 

Eliminating h and it from equations (A5) and (A6) by using equations (A10) and (All), 
equations (A5) and (A6) reduce to a set of first-order linear ordinary differential equations 
for y x and y 2 as follows: 



dyi 
dr 



where 



= { (y " 3 " r ^r) 1 + x ( JFn + c 2^2i) } yi 

+ (^12 + c 2 &aF 22 ) + C -^r (F u + c 2 &aF 2l ) - ^1 j y 2 , (A14) 

lCi<7 Ci(T 1 J 

^2 f / 2A ^2 ,x -2 2 T 1 

r— = |(e cicr + rA r ) 1 - ci<j g J^2i| yi 

+ {-M r + C/)l-x(^22 + c 2 ^ 21 )}y 2 , (A15) 

;F n = WO , (A16) 
^12 = W, (A17) 

;F 21 = 7^ WO - CoLr^ , (A18) 
F 22 = TIW, (A19) 

K = mA 1 + CoL^MoAq 1 , (A20) 

W = A (L - MxL^Mo)- 1 . (A21) 
The surface boundary conditions Ap(r = R) = are given by 

yi - y 2 = . (A22) 

The inner boundary conditions at the stellar center are the regularity conditions of the 
eigenfunctions. 



- 15 - 



REFERENCES 

Andersson, N. 1998, ApJ, 502, 708 

Andersson, N., k Kokkotas, K. D. 2001, Int. J. Mod. Phys., D 10, 381 
Beyer, H. R., & Kokkotas, K. D. 1999, MNRAS, 308, 745 
Bildsten, L., Ushomirsky, G., k Cutler, C. 1996, ApJ, 460, 827 
Bryan, G. H. 1889, Phil. Trans. R. Soc. London, A 180, 187 
Finn, L. S. 1988, MNRAS, 101, 367 

Friedman, J. L., & Lockitch, K. H. 1999, Prog. Theor. Phys. Suppl., 136, 121 

Friedman, J. L., k Lockitch, K. H. 2001, preprint (gr-qc/0102114) 

Friedman, J. L., & Morsink, S. M. 1998, ApJ, 502, 714 

Gudmundsson, E. H., Pethick, C. J., k Epstein, R. I. 1983, ApJ, 272, 286 

Kojima, Y. 1998, MNRAS, 293, 49 

Kojima, Y., k Hosonuma, M. 1999, ApJ, 520, 788 

Kojima, Y., k Hosonuma, M. 2000, Phys. Rev. D, 62, 044006 

Lee, U. 1995, A&A, 303, 515 

Lee, U., k Saio, H. 1986, MNRAS, 221, 365 

Lindblom, L. 2001, preprint (astro-ph/0101136) 

Lindblom, L., k Splinter, R. J. 1990, ApJ, 348, 198 

Lockitch, K. H., k Andersson, N. 2001, preprint (gr-qc/0 106088) 

Lockitch, K. H., Andersson, N.,& Friedman, J. L. 2001, Phys. Rev. D, 63, 024019 

McDermott, P. N., Van Horn, H. M, k Hansen, C. J. 1988, ApJ, 325, 725 

McDermott, P. N., Van Horn, H. M, k Scholl, J. F. 1983, ApJ, 268, 837 

Papaloizou, J., k Pringle, J. E. 1978, MNRAS, 182, 423 

Provost, J., Berthomieu, G., k Rocca, A. 1981, A&A, 94, 126 



-16- 

Regge, T., & Wheeler, J. A. 1957, Phys. Rev., 108, 1063 
Reisenegger, A. k Goldreich, R 1992, ApJ, 395, 240 

Richardson, M. B., Van Horn, H. M., Ratcliff, K. F., & Malone, R. C. 1982, ApJ, 255, 624 
Ruoff, J., & Kokkotas, K. D. 2001a, preprint (gr-qc/0101105) 
Ruoff, J., & Kokkotas, K. D. 2001b, preprint (gr-qc/0 106073) 
Saio, H. 1982, ApJ, 256, 717 

Schutz, B. F., & Verdaguer, E. 1983, MNRAS, 202, 881 

Thorne, K. S. 1971, in General Relativity and Cosmology, ed. R. K. Sachs (Academic Press), 
237 

Thorne, K. S. 1980, Rev. Mod. Phys., 52, 299 

Unno, W., Osaki, Y., Ando, H., Saio, H., & Shibahashi, H. 1989, Nonradial Oscillations of 
Stars, Second Edition (Tokyo: Univ. Tokyo Press) 

Yoshida, S. 2001, ApJ, 558, 263 

Yoshida, S., & Futamase, T. 2001, Phys. Rev. D, 64, 123001 
Yoshida, S., & Kojima, Y. 1997, MNRAS, 289, 117 
Yoshida, S., & Lee, U. 2000a, ApJ, 529, 997 
Yoshida, S., & Lee, U. 2000b, ApJS, 129, 353 
Yoshida, S., & Lee, U. 2001, ApJ, 546, 1121 



This preprint was prepared with the AAS IATgX macros v5.0. 



Table 1. Neutron Star Models 



Model M (M ) R (km) p c (g cm 3 ) T c (K) GM/(c 2 R) 



NS05T7 0.503 9.839 9.44 x 10 14 1.03 x 10 7 7.54 x 10~ 2 
NS05T8 0.503 9.785 9.44 x 10 14 9.76 x 10 7 7.59 x 10~ 2 
NS13T8 1.326 7.853 3.63 x 10 15 1.05 x 10 8 2.49 x 10" 1 



Table 2. Scaled Eigenfrequencies k of Fundamental r-modes 



Model 


kq(1 = m = 2) 


kq(1 — m — 3) 


2/3 x Q(0)/n 


2/3 x Q(R)/n 


NS05T7 


0.600 


0.455 


0.523 


0.645 


NS05T8 


0.601 


0.456 


0.524 


0.642 


NS13T8 


0.393 


0.309 


0.208 


0.489 
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0.1 0.2 0.3 0.4 



Fig. 1. — Scaled frequencies n = a/Q of the fundamental r-modes in the neutron star 
models NS05T7, NS05T8, and NS13T8 are plotted as functions of Cl = H/(GM/iJ 3 ) 1 / 2 for 
I' = m = 2. Note that the two r-mode frequency curves for the models NS05T7 and NS05T8 
overlap each other almost completely. 







1 1 1 1 1 1 1 


i i i i i i i 


NS05T7 


l' = m=3 
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r-mode 
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0.1 0.2 0.3 0.4 

Q 



Fig. 2. — Same as Figure 1 but for the case of m — 3. 
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r/R 



Fig. 3. — Eigenfunction i T 2 of the r-mode with V — m — 2 for the model NS05T8 at 
(l = 0.1 is given as a function of r/R. Here, normalization of the eigenfunction is chosen as 
i T 2 (R) = I. 




r/R 



Fig. 4.— Same as Figure 4 but for the model NS13T8. 
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Log(l -r/R) 

Fig. 5. — Brunt- Vaisala frequency due to the thermal stratification (solid line) and that due 
to the composition gradient (dashed line) are given as a function of log(l — r/R) for model 
NS13T8, where the frequencies are normalized by (GM/R 3 ) 1/2 . 
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0.1 0.2 0.3 0.4 
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Fig. 6. — Scaled frequencies k = a/Vt of the V = m = 2 fundamental r-modes in the neutron 
star models NS05T8 and NS13T8 are plotted as functions of Q = Q/(GM/ R 3 ) 1 / 2 . Labels 
"comp." and "temp." stand for the frequencies of the models with the buoyancy due to the 
composition gradient and due to the temperature gradient, respectively. 



-22 - 




Fig. 7. — Eigenfunctions iT 2 of the /' = m = 2 fundamental r-modes in the neutron star 
model NS13T8 at Cl — 0.2 are given as a function of r/R. Here, normalization of the eigen- 
function is chosen as % T 2 (R) = 1. The solid line and dashed line denote the r-modes in the 
models with the compositional stratification and with the thermal stratification, respectively. 




r/R 



Fig. 8. — Same as Figure 7 but for the model at & — 0.02. 



